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Models of new physics induce Kq — Kq mixing through operators having Dirac structures other 
than the "left-left" form of the Standard Model. To carry out the chiral-continuum extrapolation 
of results from numerical simulations, one needs to know the quark mass and lattice spacing 
dependence of the corresponding B-parameters in the partially quenched theory at least at next- 
to-leading order. For simulations using staggered fermions (such as that we are doing with HYP- 
smeared valence fermions on the MILC asqtad lattices) one must determine this dependence using 
staggered chiral perturbation theory (SChPT). We have calculated the required dependence in both 
SU(3) and SU(2) SChPT, working at next-to-leading order, and we give here an overview of the 
methodology and results. The SU(3) SChPT result turns out to be much simpler than that for 
the Standard Model Bk operator, due to the absence of chiral suppression for the new operators. 
The SU(2) SChPT result turns out to be closely related to that for Bk'. the chiral logarithms 
are identical, up to an operator-dependent sign. Our results are also useful for fermions with 
chiral symmetry as they provide, in the continuum limit, the partially quenched generalization of 
existing continuum results. 
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1. Introduction 

Theories of physics beyond the standard model (BSM) are often highly constrained by the ob- 
served weakness of flavor-changing neutral processes. One of the strongest constraints comes from 
— Ko mixing, both the CP-conserving and violating parts. For a given model of new physics, in- 
tegrating out heavy particles leads to a set of local AS = 2 operators with known Wilson coefficients. 
In order to determine the constraint placed on the parameters of the model, one must calculate the 
Kq — ^mq matrix elements of the AS = 2 operators. This is a task for which lattice methods are 
well suited, and several efforts are underway [JTL H|. The results will provide complementary 
information to the direct searches for new physics at the LHC. 1 

Discretization effects and unphysically large u and d quark masses are important sources of 
systematic error in these calculations, which fortunately can be largely eliminated by performing 
chiral and continuum extrapolations with effective field theory as a guide. The need for controlled, 
model-independent continuum extrapolations is particularly acute for our calculation because we 
use rooted staggered quarks, and one must extrapolate away the effects of taste-breaking and root- 
ing. To do this, one requires at least a next-to-leading order (NLO) calculation in staggered chiral 
perturbation theory (SChPT). Here we report on such a calculation, which has been recently pub- 
lished [Eh. We have presented results for both SU(3) and SU(2) SChPT, i.e. treating the strange 
quark as light and heavy, respectively. Our results are also useful for describing data obtained with 
other fermion discretizations, e.g. domain-wall fermions |[T]], for we can turn off taste breaking and 
determine the NLO predictions for a partially quenched (PQ) continuum theory in both SU(3) and 
SU(2) cases. 

Since all the technical details have been provided in Ref. we attempt here a complementary 
discussion, focusing on the essential features and eschewing most details. We also to correct some 
minor shortcomings in the original discussion, which, however, lead to no changes in the final 
results. 

BSM theories lead to AS = 2 operators with Dirac structures which are not constrained to have 
the "left-left" structure of the standard model (SM) operator. We adopt the basis 



01 = [s a y^l-y 5 )d a ][s b y^l-y 5 )d b ], (1.1) 

02 = [!"(!- 75)d a ][^(l- y 5 )d b ], (1.2) 

03 = [I"<V(1 - Y5)d a ][i , a ilv {\ - y 5 )d b ) , (1.3) 

04 = [f(\-r 5 )d a }[s b (i + r5)d h }, (i.4) 

5 = [ry )l {\-y 5 )d a }[iy )l {\ + y 5 )d b }, (1.5) 



where Gy, v = [y^,y v ]/2, and a and b are color indices. This basis was used in the highest-order 
continuum calculations of anomalous dimensions and is the most natural choice for our stag- 
gered calculation of the matrix elements, in which we match onto continuum operators using 1-loop 
perturbation theory (PT) [0]. 

Many authors, including Refs. |jl[ |[|, use an alternate basis in which O3 and O5 are replaced 

by 

3 = [f (1 " Ys)d b ] " 75)d a ] , (1.6) 

'For a recent review of this complementary information using the first lattice results see Ref. H. 
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0' 5 = [s a (l-Y5)d h }[s b (l + Y5)d a ], (1.7) 

We do not discuss the pros and cons of the two bases here, since the choice of basis turns out to 
be irrelevant for the SChPT calculation as the chiral behavior of the matrix elements of O3 and 0' 3 
(and of O5 and 0' 5 ) are identical. Thus the results presented here hold for either basis. 

In the SM, only the operator 0\ appears, its matrix element being parametrized by Bk- A large 
effort by the lattice community has led to an accurate determination of Bk- The matrix elements of 
the other operators, which we call "BSM operators", are in fact easier to calculate on the lattice. 
This is because they do not vanish in the SU(3) chiral limit, and so have simpler chiral expansions. 
Thus lattice results for the new matrix elements are following closely on the heels of those for Bk- 

As for Bk, calculating ratios of matrix elements is advantageous because many lattice sys- 
tematics and also some of the chiral logarithms cancel. Standard ratios for the BSM operators 
are 

^ = W^iTOOW ' W,^) = (5A4,-2,4/3),(l* 

where /I is the renormalization scale. Note that, unlike for Bk, the denominator does not vanish in 
the SU(3) chiral limit. Below we present and discuss our results for the chiral-continuum behavior 
of these Bj. 



2. Overview of methodology 

Our numerical calculations use a mixed-action set-up with HYP-smeared valence staggered 
quarks atop 2+1 flavors of rooted asqtad sea quarks (MILC configurations). We allow the valence 
d and s masses, labeled m x and m y , respectively, to differ from the sea-quark masses, labeled m u , 
md and m s . Although our present simulations have m u = mj, we consider the completely non- 
degenerate case in the SChPT calculation. 2 

The calculation largely follows the methodology worked out in Ref. [Jsj] for Bk in PQ SChPT 
and extended to mixed action SChPT in Ref. In most respects the calculation for the BSM 
operators is more straightforward than that for Bk because the matrix elements do not vanish in the 
SU(3) chiral limit. There is, however, an additional complication whose origin is the non-trivial 
mixing between the BSM operators in PT. 

In a continuum ChPT calculation, one would proceed by first mapping the continuum BSM 
operators into chiral operators in the effective theory, and then calculating one-loop diagrams. 
The operator mapping is done using transformation properties under the chiral group SU (Nf)L X 
SU (N/)r. To carry this over to staggered ChPT one must deal with several complications: 

• The presence of additional valence tastes implies that the lattice operators are, even in the 
continuum limit, in a partially quenched theory with an enlarged chiral group. Because 
Fierz rearrangement differs in the presence of tastes, one must introduce two down and two 
strange valence fields (D\, D2, Si and S2, where uppercase indicates a field with four tastes) 
in order to have a theory in which Wick contractions can be matched with those of the 

2 A11 masses are in physical rather than lattice units in this article. 
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desired continuum theory [Q]. One must also choose the taste of the lattice operator — with 
the standard choice being the Goldstone taste £5. 

Taste-breaking in the staggered action away from the continuum limit introduces additional 
terms in the chiral Lagrangian which are proportional to a 2 , and, in the standard power- 
counting, these terms appear at leading order (LO) along with terms proportional to m and 

The sea quarks are rooted and differ from the valence quarks. The rooting can be dealt 



with in SChPT using the replica method of Ref. [lip. The mixed action introduces further 
complications |Q] that, however, turn out to be unimportant for the BSM operators (though 
not for Bk)- 

• Taste-breaking due either to the staggered action or introduced by the use of inexact (i.e. 
one-loop) matching leads to many new chiral operators. This proliferation of operators is a 
major problem for B K [in SU(3) ChPT] but turns out to be a minor issue here as the new 
operators only contribute to analytic terms at NLO in SChPT (and not to loops). 

To deal with these complications one must follow a rather elaborate series of matching steps. The 
details are given in Ref. [||. Here we give only an overview. 

We illustrate the required steps by focusing on 0%- In step 1, we match from QCD to a PQ 
continuum theory containing two valence down quarks, d\2, and two valence strange quarks, s\^%, 
as well as their corresponding ghosts. 4 The sea-quarks are as in QCD. O2 then matches onto 

P 2 QA = 2{m - y 5 KM(l - 75)4] + - 75)d a 2 ] [4(1 - 75 Mil} • (2-D 

The two terms correspond to the two types of Wick contraction in QCD — note the [Jirf2][52^i] 
flavor structure in the second term. The overall factor of 2 compensates for an overall reduction 
by 2 in the number of Wick contractions. This matching is exact to all orders in PT. The point of 
this step is to separate the two Wick contractions, without yet having to deal with any staggered 
complications. 

In step 2 we stay in the same PQ continuum theory, but match onto operators having [jj d\] [52^2] 
form. This is in preparation for matching to the staggered theory. The result is a linear combination 
of four operators (note the different Dirac and color structures): 

of b = h m - Ys)df]m - Y5)4]+hm - y 5 )*ia - 75)4] w 

+ *3[J?(T M v(l " 75K][5 2 '<V( 1 " 75)41 +M5?Oix V (l " 75)*1<W1 - 7s)4] • (2-3) 

The coefficients bj have the form bj + abj + ... and the LO and NLO coefficents are known. 
The b^ can be obtained by a Fierz transform, while the b^p require a 1-loop calculation and one 
must make a choice of evanescent operators. 



3 This discussion extends and corrects that of Ref. g], since we have subsequently found that additional matching 
steps (steps 1 and 2 below) are required. The final results are, however, unaffected. Details will be presented in Ref. [|l3||. 
4 At this stage we can also let the valence and sea quark masses differ. 
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In step 3 we match onto the "staggered partially quenched" (SPQ) continuum theory, which 
contains sea quarks U, D and S and valence quarks D\ 2 and Si ,2 (with uppercase again indicating 
four tastes) as well as sufficient ghosts (25 in all) to cancel both the valence determinant and the 
contribution of 3 of the 4 tastes for each sea quark. Here, for simplicity, we drop the negative parity 
parts of the operator, which do not contribute to the matrix element we are ultimately interested in. 
The operator matching is exemplified by 

(2.4) 

where N t = 4 is the number of tastes. Here we are assuming external kaons of taste ^5 created with 



normalized operators. This matching is exact to all orders in PT. Each operator in the sum (2.3) is 
matched similarly and the coefficients bj are unchanged. 

We are now ready to match, in step 4, to a lattice theory with rooted staggered (asqtad) sea 
quarks % u , %d an d Xs an d (HYP) valence quarks Xdifl and Xsi2- Our continuum operator from step 
3 is in the form that is used in practice in staggered lattice calculations of matrix elements. The 
matching from the continuum PQ theory to the lattice theory is (at present) done at 1-loop in PT. 
At this order in a, a large number of lattice operators contribute to the matching, although most 
have the wrong taste and are dropped in numerical calculations. The form of the lattice operators 
that are kept is described in Ref. The errors that are introduced at this stage have the form 

lat a* spq + ^>( a )[ wrong taste ops ] _|_ 0(a z , a 2 ) [various taste ops] , (2.5) 

where = indicates that matrix elements (rather than operators) are being matched. 

In step 5 we match onto the Symanzik effective continuum theory describing the long distance 
degrees of freedom of the lattice theory. Lattice symmetries are used to determine corrections to 
the action and operators, and factors of a 2 are now explicit. 

Finally, in step 6 we match onto the chiral effective theory (SChPT), using the chiral transfor- 
mation properties of the action and operators. This process is straightforward, though tedious. A 
crucial step in this matching (which we do not have space to enter into) is that certain symmetry 
properties of the original operator percolate through the multiple steps and constrain the low-energy 
constants (LECs) that multiply the chiral operators, reducing the number of independent constants. 
The details of this are slightly modified from the discussion in Ref. ^ by the addition of steps 1 
and 2, but the final conclusion is the same. 

With the chiral form of the operators in hand, the requisite one-loop calculation is now straight- 
forward. We have done this in SU(3) SChPT, and obtained the SU(2) result using the method 
explained and justified in Refs. [^, ||]. 

3. Results and conclusions 

The general form of the NLO result in SChPT is 



B i= B LO 



1 + SBf al + 8B l ° g 



(3.1) 



The analytic terms involve five unknown LECs at NLO: 

<5fi™ al = C ji (m x + m y ) + Cj 2 (m u + m d + m s ) + cpa 2 + c^a 2 + Cj 5 a 2 . (3.2) 
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Here we have shown the SU(3) result; for SU(2) the m y and m s terms are absent (and the coefficients 
are different). The first two terms are present in the continuum, while the remaining three represent 
lattice artifacts. The a 2 term comes from standard discretization errors (which are not suppressed 
by a since our valence fermions are not fully improved). The a 2 a = a 2 a 2 term comes from taste- 
breaking discretization errors, which are numerically enhanced such that they are comparable to 
standard &(a 2 ) errors. Finally, the a 2 term arises because we use 1-loop operator matching — it 
would be absent if we used non-perturbative renormalization. 

These analytic terms are similar to those for B K and one can deal with them in a similar fashion. 
In a chiral extrapolation at fixed a the last three terms are constants and can be absorbed into B 1 - . 
These terms then enter into a (fairly complicated) continuum extrapolation. 

The non-analytic terms are completely predicted in terms of the pion decay constant in the 
appropriate chiral limit or fi) and the masses of the valence and sea pions. They lead to 
curvature in the predicted chiral behavior. We show explicitly only the SU(2) result for Nf = 2 + 1 
(which is the most useful in practice): 



8B l ° s = ±— i - 

3 (47T/ 2 



1 + \ {mi-i + {Mlj-M 2 xx .j)i{M 2 xx .j)} 



B 



(3.3) 



Here M xx - B is the LO mass of the valence xx pion of taste B and M K -j is the LO mass of the taste- 
singlet sea-quark pion. The chiral logarithms are £(X) =X\n(X/iJ. 2 ) and £(X) = —d£(X)/dX. The 
+ sign holds for j = 2 and 3, while the — sign applies for j = 4 and 5. 
Key features of this result are: 

• The chiral logarithms for j = 2 and 3 (and for j = 4 and 5) are identical because the 
chiral transformation properties of these pairs of operators are the same Jl2[] . This holds 
generally — i.e. also for SU(3) SChPT, and irrespective of degeneracies in sea-quark masses. 

• The j = 2,3 and j = 4, 5 cases differ only in sign. This relationship does not hold for the 
SU(3) result. 

• The chiral log for j = 2,3 is identical to that for Bk ■ This identity does not hold for the SU(3) 
result. 

• The result is unaffected by the use of a mixed action — this only impacts some coefficients 
in 5By nal as well as the values of the pion masses and decay constant. This conclusion holds 
also for the SU(3) result (but not for Bk)- 

In Ref. ^ we present also the SU(2) result for non-degenerate sea quarks, and the SU(3) result 
both for Nf = 1 + 1 + 1 and Nf = 2+ 1. We only note here that the SU(3) chiral logarithms for the 
BSM operators are much simpler than those for Bk- The latter involve 13 additional LECs, whereas 
the former involve none. This simplicity can be traced back to the fact that the matrix elements of 
the BSM operators do not vanish in the SU(3) chiral limit. 

Some time ago, Becirevic and Villadoro (BV) calculated the chiral logarithms for the BSM 



operators in continuum, unquenched SU(3) ChPT [12] (results which we use as a check on ours 
in the appropriate limit). They noted that the chiral logarithmss cancel in the ratios B^jB^ and 
B4/B5. They dubbed these ratios "golden" and advocated their use in minimizing errors in chiral 
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extrapolations. Our results allow us to extend the set of golden combinations. First (and this is 
almost trivial) the two golden ratios of BV remain golden in the PQ theory, and in SU(2) ChPT. 
Second, the ratio B2/BK becomes golden in SU(2) ChPT. Finally, the product B2B4 is also golden 
in SU(2) ChPT. Neither of these new combinations is golden in SU(3) ChPT, although there is a 
partial cancellation of chiral logarithms (and so in BV's parlance they are "silver" combinations). 

These golden combinations could be particularly important for a calculation using staggered 
fermions. This potential importance arises because chiral logarithms are responsible for a signifi- 
cant fraction of the taste breaking in the Bj (since the pion masses entering the logarithms are taste 
dependent). Thus the chiral behavior of the golden combinations may be smoother and easier to 
fit. One need only fit one of the five B-parameters including chiral logarithms, determining the 
other four using golden ratios or products. We are presently testing this approach in our companion 
numerical calculation 

Acknowledgments 

W. Lee is supported by the Creative Research Initiatives Program (2012-0000241) of the NRF 
grant funded by the Korean government (MEST), and acknowledges support from KISTI super- 
computing center through the strategic support program for the supercomputing application re- 
search [No. KSC-201 1-G2-06]. S. Sharpe is supported in part by the US DOE grant no. DE-FG02- 
96ER40956. 

References 

[1] P. A. Boyle, N. Garron and R. J. Hudspith, Neutral kaon mixing beyond the standard model with 
re/ = 2 + 1 chiral fermions, Phys. Rev. D 86, 054028 (2012) [arXiv: 1206 . 5737 [hep-lat] ]. 

[2] V. Bertone et al., Kaon Mixing Beyond the SM from Nf=2 tmQCD and model independent constraints 
from the UTA, arXiv: 1207 . 1287 [hep-lat] . 

[3] J. A. Bailey et al., SWME Collaboration, these proceedings: Beyond the Standard Model corrections 
to K° -K° mixing, arXiv: 1211. 1101 [hep-lat] . 

[4] F. Mescia and J. Virto, Natural SUSY and Kaon Mixing in view of recent results from Lattice QCD, 

arXiv: 1208 . 0534 [hep-ph]. 

[5] J. A. Bailey, H. -J. Kim, W. Lee and S. R. Sharpe, Kaon mixing matrix elements from 

beyond-the-Standard-Model operators in staggered chiral perturbation theory Phys. Rev. D 85, 
074507(2012) [arXiv: 1202 . 1570 [hep-lat]]. 

[6] A. J. Buras, M. Misiak and J. Urban, Two loop QCD anomalous dimensions of flavor changing four 
quark operators within and beyond the standard model, Nucl. Phys. B 586, 397 (2000) 

[hep-ph/0005183]. 

[7] J. Kim, W. Lee and S. R. Sharpe, One-loop matching of improved four-fermion staggered operators 
with an improved gluon action, Phys. Rev. D 83, 094503 (2011) [arXiv: 1102.1774 
[hep-lat ] ]. 

[8] R. S. Van de Water and S. R. Sharpe, Bk in staggered chiral perturbation theory, Phys. Rev. D 73, 
014003 (2006) [hep-lat/0507012]. 



7 



Chiral extrapolation of matrix elements ofBSM kaon operators 



Stephen R. Sharpe 



[9] T. Bae, et al, Bk using HYP-smeared staggered fermions in Nf = 2 + 1 unquenched QCD, Phys. Rev. 
D82, 114509(2010)[arXiv:100 8.517 9 [hep-lat]]. 

[10] W. -J. Lee and S. R. Sharpe, Partial flavor symmetry restoration for chiral staggered fermions, Phys. 
Rev. D 60, 114503 (1999) [hep-lat/ 9 905023]. 

[11] C. Aubin and C. Bernard, Pion and kaon masses in staggered chiral perturbation theory, Phys. Rev. D 
68, 034014(2003) [hep-lat /0304 014]. 

[12] D. Becirevic and G. Villadoro, Remarks on the hadronic matrix elements relevant to the SUSY 
K -K^ mixing amplitude, Phys. Rev. D 70, 094036 (2004) [hep-lat / 4 8 02 9]. 

[13] J.A. Bailey, H.-J. Kim, W. Lee and S.R. Sharpe, in progress. 



8 



